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We investigate many-body localization in the presence of a single-particle mobility edge. By 
considering an interacting deterministic model with an incommensurate potential in one dimension 
we find that the single-particle mobility edge in the noninteracting system leads to a many-body 
mobility edge in the corresponding interacting system for certain parameter regimes. Using exact di- 
agonalization, we probe the mobility edge via energy resolved entanglement entropy (EE) and study 
the energy resolved applicability (or failure) of the eigenstate thermalization hypothesis (ETH). Our 
numerical results indicate that the transition separating area and volume law scaling of the EE does 
not coincide with the non-thermal to thermal transition. Consequently, there exists an extended 
non-ergodic phase for an intermediate energy window where the many-body eigenstates violate the 
ETH while manifesting volume law EE scaling. We also establish that the model possesses an infi¬ 
nite temperature many-body localization transition despite the existence of a single-particle mobility 
edge. We propose a practical scheme to test our predictions in atomic optical lattice experiments 
which can directly probe the effects of the mobility edge. 


Thermalization, a commonplace phenomenon in vari¬ 
ous physical settings, can naturally fail in isolated dis¬ 
ordered quantum interacting systems, making standard 
concepts of quantum statistical mechanics invalid. The 
fundamental theoretical underpinning of thermalization 
in quantum systems has been postulated in the form of 
the eigenstate thermalization hypothesis (ETH) [U [2]. 
Recently, it has been shown using perturbative argu¬ 
ments that the presence of interaction and disorder in 
a closed quantum system could lead to many-body lo¬ 
calization (MBL) [3] with such an interacting quantum 
MBL state being non-thermal. 

A hallmark of MBL is its violation of the ETH [2], 
where a local subsystem fails to thermalize with its en¬ 
vironment [4]. MBL has now been established non- 
perturbatively in lattice models with hnite energy den¬ 
sity, where numerical evidence points towards the exis¬ 
tence of MBL all the way to infinite temperature [SJE]. 
Further numerical work [THS] and a rigorous mathemat¬ 
ical proof [To] for the existence of the MBL phase have 
mounted compelling evidence for the existence of such a 
‘finite-temperature’ MBL phase which eventually gives 
way to an extended phase at strong enough interac¬ 
tion. Although much of the MBL work has focused on 
the interacting one dimensional (Id) fermionic Anderson 
model with random disorder m (and closely related spin 
models), it turns out that MBL also exists without any 
disorder li m US] for the Aubry-Andre-Azbel-Harper 
(AAAH) model [T^16| . which is a non-random Id model 
with a quasiperiodic onsite potential. We emphasize that 
neither Id Anderson model nor AAAH model manifests 
a single-particle mobility edge (SPME). 


introducing a SPME allows one to study how localized 
and delocalized eigenstates will interact, thus introducing 
qualitatively new physics. There are several determinis¬ 
tic Id incommensurate models with SPMEs in the liter¬ 
ature HM], which can be adapted for studying MBL 
in the presence of a SPME. 

We consider a recent generalization |23| of the Id 
AAAH model with an analytical expression for the 
SPME, which enables us to study the interplay of many- 
body effects and the SPME in a controlled fashion. Since 
the MBL phase is a property of all eigenstates, the pres¬ 
ence of a mobility edge adds a new dimension to the 
problem as both localized and delocalized single-particle 
orbitals are now present in the problem. Using exact di- 
agonalization we find for certain parameter regimes of the 
model: (1) The existence of a many-body mobility edge 
(El) characterized by the area to volume law scaling of 
entanglement entropy (EE). (2) A distinct energy scale 
(Et) that separates a thermal (i.e. ergodic) and non- 
thermal region in energy, which is established by directly 
considering an ETH violation based on the criterion in 
Ref. [1112121] • (3) Our results suggest Et ^ El and con¬ 
sequently the existence of a non-ergodic regime with vol¬ 
ume law EE scaling between Ep and El- All three of our 
findings are completely novel differing drastically from 
previous studies suggesting a sharp many-body mobility 
edge piiaiis]. To guide future experiments that could 
probe our predicted mobility-edge physics, we present a 
realistic scheme with a straightforward modification to 
the existing experimental setup P7H35] . 

The model we consider is a generalized Aubry-Andre 

(GAA) model [23], = iLo + 


In the absence of a SPME, interactions act on the Fock 
space of Slater determinants of either completely local¬ 
ized or delocalized single-particle eigenstates. Therefore, 
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Figure 1. Energy dependent ergodicity for interacting fermions, (a), A scenario for many-body mobility edge physics as we 
have found for the GAA model, (b) The bipartite Renyi entropy. The eigenstates above El have extensive EE, and are thus 
extended; whereas the states below El exhibit EE of area law scaling and are localized, (c) The energy dependence of an 
observable 0{E) (Eq. §) with L = 30. The inset shows the fluctuation var[Cl] (see main text) for different system sizes, using 
the same plot scheme as in (b). The fluctuation of O for states E > Et is small and gets significantly larger for states E < Et- 
The system is expected to be thermal (non-thermal) above (below) Et- (The slight increase of var[Cl] near E/N ~ 1.2 is an 
artifact from the small number of states close to the spectra-edge.) Our numerical results suggest El < Et- In this plot, the 
filling is fixed at 1/6, and we use X/t — 0.3, V/t = 1, and a = —0.80, and we average over </ for better statistics [30]. In the 
calculation for system size L — 30, we use inverse Lanczos and target 500 interior eigenstates. 


where Cj is a fermionic annihilation operator, rij = CjCj, 
and the tunneling t is the energy unit throughout. We 
focus only on the fermionic case here. We consider a G 
(—I, I), with the AAAH model corresponding to a = 0. 
In the non-interacting limit, the GAA model with an 
irrational wavenumber q (we fix q = 2/(1 -|- -s/S), with no 
loss of generality), has a SPME |23j at ae = 2sgn(A)(|t| — 
|A|) jSni- In this model, the particle number ^ 

is conserved. 

Interaction effects on localization and thermalization 
in the presence of SPME- — We study interaction ef¬ 
fects on localization and thermalization for the model 
Hamiltonian H using exact diagonalization. In the 
AAAH model (a = 0), the non-interacting many-body 
wavefunction is a Slater determinant of all localized or 
all extended single-particle orbitals. This results in the 
interacting AAAH model having all many-body states 
either localized and non-thermal or extended and ther¬ 
mal |8]. 

However, for the non-interacting GAA model (with 
a ytz Of there are more possibilities originating from the 
SPME, where the Slater determinant can be composed 
of both localized and extended single-particle orbitals. 
Adding interactions to such a system may result in richer 
many-body states where the localization and thermaliza¬ 
tion properties may be qualitatively different from the 
a = 0 case. To this end, we employ separate diagnos¬ 
tics to study localization and ergodicity without making 
the common assumption that thermalization and delocal¬ 
ization must necessarily be intrinsically connected in an 
interacting system. To investigate the localization prop¬ 
erties, we cut the lattice at site I, which divides it into 
two subsystems A and B, we then calculate the energy re¬ 
solved Renyi entropy 82 ( 1 ) = — log(Trp^) of A involving 
lattice sites 1, 2, ..., I, whose reduced density matrix is 
obtained by tracing out region B at the other sites {l + l, 


1 + 2,-. -, L) |3T]. The EE scaling reliably tracks the local¬ 
ization transition, where localized and delocalized many- 
body states are quantified by the area law (S '2 ^ 
and volume law scaling {S 2 ^ L'^) respectively |4l|8|. To 
understand the thermalization features we calculate the 
observable 0{E), 

= ( 2 ) 

with IdlB) a many-body eigenstate. The large fluctuation 
in 0{E) among eigenstates that are nearby in energy is 
a signature for the violation of the ETH 

We begin by focusing on the MBL transition as a func¬ 
tion of energy for fixed model parameters. In Fig. [^we 
show the energy resolved 82(1 = L/ 2 ) for various system 
sizes. We find eigenstates with an energy below a certain 
value El are localized with an EE that obeys area law 
scaling (5'2(T/2) ^ L°) |3], whereas the eigenstates with 
an energy above El exhibit a volume law scaling of the 
EE {S 2 {L/ 2 ) ~ L) and are thus extended. We define 
El where S 2 {L/ 2 ) splays out in system size as shown in 
Fig-^- Thus El defines the many-body mobility edge, 
which separates states with an area law scaling from ex¬ 
tended states with volume law scaling. Although the ex¬ 
istence of the many-body mobility edge El Ai our model 
is already a significant result, below we discuss the key 
issue of whether El also defines the ergodic properties of 
the interacting system. 

We now come to the thermalization properties, which 
are captured by the energy resolved observable 0{E) as 
shown in Fig. [^. The fluctuations of 0{E) within a 
narrow energy window are quantified by their variance, 
denoted as var[(!I] [30]. As shown in Fig. there is a 
clear energy threshold Et, that separates two qualita¬ 
tively different regimes. In the energy window {E > Et), 
the fluctuations of O among nearby eigenstates are small, 
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for which the ETH is satisfied and the eigenstates are 
thermal [24] . The spread of the observable broadens out 
for energies E < Et-, where the fluctuations are signif¬ 
icantly larger, leading to a violation of ETH [53] . We 
emphasize here that the thermal to non-thermal transi¬ 
tion is unique to interacting systems, and is absolutely 
absent without interactions [531 132| . 

Our numerical results suggest that the MBL and ther¬ 
mal transitions in energy do not coincide, El ^ Et 
(see Fig. [^. In particular, our numerics imply that 
Et > El, and as a result there exists an energy win¬ 
dow {El < E < Et) where the many-body states are 
non-ergodic (violate the ETH) but remain extended (vol¬ 
ume law scaling of EE). Therefore, we conclude that in 
this interacting many-body system with SPME two criti¬ 
cal energy scales El and Et exist, which is qualitatively 
different from the scenario of a sharp many-body mo¬ 
bility edge [5S]. We emphasize that our results are com¬ 
pletely distinct from the model in the absence of a SPME 
(a = 0), which has no many-body mobility edge and all 
of the eigenstates are either thermal and delocalized or 
non-thermal and localized |30j . Thus our numerical re¬ 
sults point to the existence of a non-ergodic extended 
regime defined as El < E < Et- The possible exis¬ 
tence of a non-ergodic extended (i.e. metallic) phase in 
the vicinity of the MBL transition has been speculated 
in the literature with no concrete examples |33H35j . 

Non-interacting many-body states. — In the follow¬ 
ing, we develop a physical intuition for the ob¬ 
served many-body mobility edge. The non-interacting 
many-body states are trivially non-ergodic and violate 
ETH [35], as shown in Fig. [^c) for E = 0 where the 
energy dependence of 0{E) manifests large fluctuations 
among eigenstates that are nearby in energy [53]. We 
will discuss the non-interacting limit to gain insight into 
the emergence of the non-ergodic extended phase in the 
interacting case. 



Figure 2. Localization and non-thermalization of non¬ 
interacting fermions with a SPME. Here we simulate 5 par¬ 
ticles in 30 sites. We take the GAA model with A/t = 0.3, 
a = —0.8, with mobility edge £o/t = —1.75. (a) illustrates 
different possibilities of many-body states. As shown the 
partially-extended state as marked by 3 could have approxi¬ 
mately the same energy as the localized one marked by 2. (b) 
shows the entanglement scaling for the three types of states, 
localized, extended and partially-extended. The partially- 
extended states exhibit extensive EE, similar to the extended 
ones. In (c) we show the energy dependence of O (Eq. ©)• 


Without interactions, the many-body eigenstate of N 
fermions is a product state of N single-particle orbitals. 
In the presence of SPME, there are three different ways of 
constructing many-body states (Fig. [^a)): (i) all parti¬ 
cles put in localized single-particle orbitals, (ii) all parti¬ 
cles in extended orbitals, and (iii) some particles in local¬ 
ized and others in extended orbitals, which respectively 
give localized, extended and partially-extended many- 
body states. The partially-extended states have exten¬ 
sive EE (Fig. ib)). (We mention that the partially- 
extended states would appear localized [30] from the per¬ 
spective of the normalized participation ratio [S].) Such 
partially-extended many-body states lead to important 
consequences. Consider a model with single-particle en¬ 
ergies ei < €2 < ■ ■ ■ < Cl having a SPME e^*, such that 
the states with em<m* (£m>m*) are localized (extended). 
The lowest energy for a fermionic many-body partially- 
extended state is Ea = em*+i + X]m=i "bhe highest 
energy of a localized state is Eb = Y.Zi=m*-N+i ^rn- For 
a general Hamiltonian, Eb > Ea is the most typical sce¬ 
nario [37]. The many-body states with energy <(>) Ea 
are completely localized (extended, partially-extended or 
localized); the states above Eb are extended or partially- 
extended. In the energy regime Eb > E > Ea, however, 
localized and partially-extended states are coexisting by 
virtue of the SPME in the spectrum which enables the 
existence of this mixed intermediate energy regime. 

Putting the non-interacting and interacting results all 
together, a physical picture naturally emerges. Inter¬ 
action effects on the extended many-body states lead 
to thermal behavior, whereas for localized states, inter¬ 
actions make them non-thermal. The most interesting 
case is the coexistence of partially-extended and localized 
states, where interactions could stabilize a non-thermal 
extended phase. Such a physical scenario is possible and 
seems to be consistent with our numerics. 

The MBL phase in the GAA model. — Past studies 
have established the existence of a MBL phase (with all 
eigenstates localized) in the AAAH model [3[ [35] at in¬ 
finite temperature. The infinite temperature limit is de¬ 
fined by averaging observables over all energy eigenstates 
(i.e. with a thermal weighting factor equal to unity). For 
the non-interacting GAA model with a SPME, the aver¬ 
age is performed over localized, extended, and partially- 
extended states, and thus leading to extended behavior, 
e.g., the averaged EE obeys volume law. For the inter¬ 
acting case, in contrast to MBL studies in the absence 
of SPME where interactions make the system more ex¬ 
tended, we find that interactions could actually stabilize 
the infinite temperature MBL phase in a considerable 
parameter region of the model in Eq.Q. 

To capture the MBL transition at infinite temperature, 
in addition to the EE scaling (averaged over all eigen¬ 
states), we present the level statistics, an established di¬ 
agnostic for MBL [51 [S]- We consider the dimensionless 
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Figure 3. MBL phase at half filling and infinite temperature, 
(a), (b) and (c) show the Renyi entropy versus the incommen¬ 
surate potential strength, (d), (e) and (f) show the averaged 
adjacent gap ratio. In (a) and (d), we show the interaction 
dependence with a = —0.8. The numerical results indicate 
the interacting system is in the MBL phase (all states local¬ 
ized) when A is larger than a certain critical value (Ac ~ 0.8). 
In particular (a) explicitly shows interaction effects make the 
GAA model more localized, (b) and (e) show the a depen¬ 
dence with V/t = 1. As we increase a, Ac becomes smaller. In 
(c) and (f), we show the system size dependence with V/t = 1 
and a = —0.8. 

adjacent gap ratio, 

Tn = min(d„,(5„+i)/max((5n,i5n+i), 

where = i?„+i — En- We calculate its average r = 
^ Hilbert space dimension), to locate 

the MBL phase [S]. 

As shown in Fig. when the incommensurate poten¬ 
tial strength A is weak, the EE is extensive, a signature 
for the system being extended. The average adjacent gap 
ratio is f « 0.53, which implies that the energy spectra 
satisfy Gaussian orthogonal ensemble level statistics, and 
the many-body phase is delocalized EliaEH]- When A is 
above a certain threshold Ac(q!, V), the model undergoes 
the MBL transition. In this parameter regime, the EE 
obeys area law scaling (Fig. and the gap ratio be¬ 
comes f = 0.39 (Fig.|^), which is consistent with eigen¬ 
states satisfying a Poisson distribution, and the model 
is in the MBL phase [2 [g 0 EH]. We emphasize that 
for A > Ac in the presence of delocalized single-particle 
orbitals, the MBL phase can still be stabilized by inter¬ 
actions. For example, with \/t = 1.5 and a = —0.80, 
although the non-interacting case is not completely lo¬ 
calized due to the SPME, the interacting system is com¬ 
pletely localized as implied by the EE (Eig. adjacent 
gap ratio (Eig. Hi). and their system size dependence 
(Fig-ic ,f)). In addition, starting from a = 0 (with no 
SPME) as a —>■ —0.99, where more localized orbitals are 


mixed in, we find both S 2 {L/ 2 ) and f decrease as dis¬ 
played in Fig. [^b,e). To conclude, we have established 
the existence of an MBL phase at infinite temperature in 
the presence of SPME. 

Experiment. To study the MBL phase in the AAAH 
model, a two-component Fermi gas of atoms has 
been recently confined in a Id superlattice with optical 
potential, Vocos‘^{kx) + Vi cos‘^{k'x) [ 37 ], with k' incom¬ 
mensurate to A:. To investigate our predicted mobility 
edge physics, we propose to add an additional poten¬ 
tial, V 2 cos^(2A:'a;). Choosing Vg = bE^, Vi = O.lSEr, 
and V 2 = 0.026Er, {Er is the single-photon recoil en- 
ergy) |39j . the non-interacting model in Eq. ([B with 
X/t = 1, and a = 0.2 is approximately realized j.'fflj . with 
localization properties described in Ref. |23|. Preparing 
an initial state with its average energy in the non-thermal 
extended region, its unitary evolution would provide di¬ 
rect observation of nontrivial relaxation of an interacting 
many-body state in presence of SPME jM] [JD]. As the 
numerical simulations on classical computers are limited 
in terms of system size, the quantum simulator, atoms 
in the optical lattice, would clarify the thermodynamic 
limit of our proposed non-ergodic and MBL phenomena. 

Conclusions:- In summary, we have shown the single¬ 
particle mobility edge and interactions result in a many- 
body mobility edge. A central new result here is the 
existence of two characteristic many-body energies (i.e. 
El and Et) in general in a system with a correspond¬ 
ing SPME, which separate localized and extended states 
{El) and non-ergodic and thermal states (Et). Our 
numerical results (within our numerical accuracy and 
within the finite size limitations) suggest El < Et, 
which allows for the possibility of non-ergodic delocal¬ 
ized many-body states (i.e. a non-ergodic metal) as a 
strange new intermediate phase of quantum matter. We 
expect our findings to generically apply to systems with 
SPME, specifically, to the three dimensional interacting 
disordered Anderson model. 

Note added. After completing our work we be¬ 
came aware of a complementary and independent recent 
study |41j of many-body localization in systems with mo¬ 
bility edges. 
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Supplementary Material: 

Many body localization and quantum non-ergodicity in a model with a 

single-particle mobility edge 


Fraction of delocalized states as a function of {a, A). 


In this section we compute the fraction of localized states as a function of A for different values of a considered in 
the main text. In Supplementary Figure]^ we plot the ratio of the delocalized states to the total number of states 
as a function of X/t for different values of a= (0, -0.2, -0.8, -0.99). For a = 0, all the states are either delocalized or 
localized as a function of X/t (with X/t = 1.0 being the critical point). For a ^ 0, there is an energy dependence to 
the fraction of delocalized states and the many step structure is in part due to the gap in the single particle energy 
spectrum. Thus the parameter a allows us direct control of the fraction of single particle localized (delocalized) states 
participating the the non-interacting many body state given by the slater determinant. 



Supplementary Figure 4. Number of delocalized state to the total number of states as a function of X/t for different values 
of a= (0, -0.2, -0.8, -0.99) 


Details of localized states with interactions 

In Supplementary Figure we provide more results about the observable 0{E) = showing its behavior 

in completely localized and extended phases, where the non-ergodicity is consistent with the entanglement entropy 
scaling. 

In numerics, the fluctuation of observables {0{E)) between nearby eigenstates, is quantified by the standard devi¬ 
ation of 0{E), var[0], over some small finite energy range [E — 6E/2,E + SE/2] (we choose SE/{Nt) = 0.02 in this 
work). 













Supplementary Figure 5. Completely localized and extended phases in the AAAH model, (a) and (b) show the fluctuation 
in the observable 0{E) in the extended and localized phases, respectively, (c) and (d) show the corresponding Renyi entropy 
scaling. Here, in calculating entanglement entropy, we average over all eigenstates. In this plot we fix V/t = 1, a = 0, N = 8, 
L — 16, and we use X/t = 0.3 and X/t = 1.5, respectively for the extended and localized phases. In the localized phase (b,d), 
we have large fluctuation in 0{E) for nearby eigenstates, a signature for non-ergodicity, and the entanglement entropy obeys 
area law. In the extended phase (a,c), the fluctuation in 0{E) is suppressed, and the entanglement entropy obeys volume law. 


In averaging over the phase cj) (see Eq.(l) in the main text), we choose 10 realizations with c^/tt = {0,0.2,... 1.8}. 
We first calculate var[(!I] for each (^-realization and then take the average. As the energy spectra for different (f) are 
not identical, we actually average over different data points within an energy window set by SE'/{Nt) = 0.04. 


Detailed description about experimental realization 


In experiments to probe the MBL phase in the AAAH quasiperiodic model, a Id superlattice is engineered with 
counterpropagating laser beams, and the resultant optical potential reads VoCOS^(A:a;) + ViCOS^{k'x) [27], with k' 
incommensurate to k. With Vb 3> Er {E^ is the single photon recoil energy), the tight binding model is a valid 
description for the experimental system. Further requiring 0 < \Vi\ ^ Vb, the AAAH model is reached. Taking 
Vb = 5Ej., the parameters in the AAAH model are then given by {t « 0.065Er, X/t « iVi/t). 

Considering the limit |a| ^ 1 of the GAA model (see Eq. (1) of the main text), the on-site potential becomes 
~ 2A(cos(27rjA:-|-((>)-|-Q!COS^(27r}fc-l-^)). It follows that to implement the GAA model we need an additional potential, 
Vb cos^(2fc'a;). Using Vb <C Vi, the GAA model is approximately realized with a estimated to be a « Vb/Ui. Choosing 
X = t, the small a limit supports both of sharply localized and delocalized wavefunctions (see Ref. [23]), and there 
is a well-defined single-particle mobility edge. It is worth noting here that the a <C 1 limit may or may not be an 
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experimental sweet spot to probe the mobility edge physics, as it requires fine tuning of the laser strength to get 
X = t. Searching for the parameter region that is most suitable for experiments is not the focus of this work and is 
thus left for future. 


Normalized participation ratio for non-interacting many body states 


To describe a many body state, we can choose a Fock state basis 

with N the particle number. We will take an ordering ji < j 2 < ■ ■ ■ Jn- The many body wavefunction in this basis is 


1^) = X! (3) 

jl<j2<--<jN 

from which the NPR is defined to be [8] 

^(1^)) = -[4-[4^, (4) 

^h<h< - <0N 

with Vh the dimension of the many body Hilbert space, which is ^ ^ 

Without interaction, the many body eigenstate is simply a product state of single-particle orbitals, namely 

. •. V’mjv l^ac). (5) 

First, as constructed in Eq. ([^, the many body wavefunction is given as 




free 

{iij2,...ijv} 


= ^sgn(P)V'„,0p[i])... 

'4^mN {jp[N]), 


( 6 ) 


where P labels all permutations. Now we get 

/4(|mi,m2,...,m4)) = ^ { 31 , 32 ,. 

31<32<..<3N 

= Sgn{PaPbPcPd) Y [V'mi(jP„[l])---V^m«(jP„[Afl)] 

Pa,Pb,Pa,Pd 31<32<...<3N 

X [V'miOpdl])---^m«(jPdAf])] [V'miOpqi])---V’m«OPc[Af])] [V'mi (jP^ [ 1 ] ) • • ■ V'mw [AT] )] 

= Y Sgn(T’aA'^c) Y [V'miOp'[l])---^m«(jPHAf])] Op'[1]) • ■ ■ Op'[A f]) 

PLPhPa 3l,32,...,3N 

X [V'mi 0 p^[ 1 ]) • ■ • V'miv (jP.^[Af])] [V'mi (jl) ■ • • Un)] ■ ( 7 ) 

If one orbital, say mi, is localized, we take '4’miij) ~ ^j.i (ap to a local unitary transformation), and we get 

h= Y SgniPaPbPc) Y [^^2 Upa [ 2 ] ) • ■ • 0 p„ [N] )] [^*m2 Pb [ 2 ] ) ' ' ' Up, [AT] ) 

Pa,Pb,Pa 32,33 ■■ .,3 N 

X '>Pm2 U P,, [ 2 ] ) • ■ • ’^rriK U P,. [AT] )J ['07712 (j 2 ) ■ ■ ■' 0 mjv (jTv)] I (^) 

where P means a permutation among the numbers (2,3,4,..., N). It follows that we have a particle number reduction 
relation 

/4(|TOi,m2, ... ,mjv)) = / 4 (|to 2 ,P 13 , ■ ■ ■ ,TOAr))- (9) 

The Hilbert space dimension for N — 1 fermions on L sites is ^ know that at low filling 

/4(|m2,TO3, . . . ,mAr)) > 


L 

N-1 


-1 


> 


-1 
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Then we have rj{\mi, m 2 ,..., itin)) —>■ 0, as L —> 00 , if we fix the particle number. 

Suppose we have a number, iVioc, of single particle localized orbitals, mi, m 2 , ..., and rriNi^^, then it follows from 
the reduction relation (Eq. @) that 

/4(|mi,m2,... ,mAr)) =/4(|mAri„^+i,... mAr)). (10) 

Then the NPR value satisfies 

f N \ 

7?(|mi,m2,...,mAf)) < K j ■ (H) 


/ \ Nlac 

Considering the thermodynamic limit, L —> 00 while keeping N/L = ptot, V < ( • To have p —>■ 0, 

the nnmber of localized single-particle orbitals composing the many body states is required to be also extensive, i.e., 
Moc = PiocL (in principle, pioc could be an infinitesimal number), as we take the thermodynamic limit. Our conclusion 
is, if the many body state is composed of some finite (extensive) fraction of localized single-particle orbitals, then its 
NPR value vanishes. The reverse is also true at low filling. The most general case shall be discussed elsewhere. 

The above analysis based on NPR shows that a many body state composed of both localized and extended orbitals 
behaves like a localized state, which is opposite to using the entanglement entropy scaling criterion. Nonetheless, 
when the system is in a MBL phase with all states localized, NPR and entanglement entropy give consistent results 
(see Supplementary Figure]^ and Figure 3 in the main text). For example, with \/t = 1.5, a = —0.8, and V/t = 1, 
the three different diagnostics, including NPR, entanglement entropy and adjacent gap ratio consistently support all 
many body states are localized. 





Supplementary Figure 6. The normalized participation ratio as a function of A for the GAA model, (a) shows the interaction 
dependence with a = —0.8. (b) shows the a dependence with V/t = 1. (c) shows the finite size dependence with V/t = 1 and 
a = —0.8. We choose /> = 0 and N/L = 1/2 in this plot. 


Similarity of the GAA to 3d Anderson model 

In Supplementary Figurej^ we show the Von Neumann entropy for single-particle states in Id GAA and 3d Anderson 
models. The behavior of localized and extended states appears similar in these two models. We suspect that the 
many body localization phenomena in the Id GAA model could resemble that in the 3d Anderson model. Detailed 
comparison of these two models will be discussed elsewhere (work in progress). 
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Supplementary Figure 7. Von Neumann entropy for single-particle extended and localized states, (a) and (b) show the 
results obtained from the Id GAA and 3d Anderson models m, respectively. For the GAA model, we choose X/t = 0.3 and 
a = —0.8, and the phase (j) is averaged over. For the 3d Anderson model, we choose the disorder strength to be sixteen times of 
the nearest neighbor tunneling, and we average over two thousand disorder realizations. For both localized and extended states, 
we simply pick certain energies, and the qualitative behavior of the entanglement entropy does not depend on the details. 





